clc;
clear;
syms x;
a=sym('a%d',[1,5]);
A=hankel(a);
result=collect(charpoly(A,x),x)
latex(result)